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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2647 

IMPLICATION OF THE TRANSPORT EQUATION FOR THE 
SEMTRMPIRICAL TREATMENT OF SHIELDS 
By Philip Schwed 


SUMMARY 

The semiempirical method of dealing with shields hy treating them 
as composed of layers is revised and extended hy taking the angular dis- 
tribution of the radiation into account and by making use of the trans- 
port equation. It is shown that breaking up the ranges of direction 
and energy of radiation into finite intervals is appropriate where mul- 
tiple scattering with angular deviation does not dominate, and a pro- 
cedure is given for using data referring to very thin layers to calculate 
the parameters for finite layers in this case. In the contrary case 
where multiple elastic scattering is dominant, it is shown that the num- 
ber of unknowns to be considered is decreased if the angular distribution 
is represented by an appropriate expansion in terms of Legendre poly- 
nomials and the Laplace transformation of the coefficients is taken with 
respect to energy. The method is illustrated by consideration of the 
effect upon neutrons of a thick nonabsorbing shield of high atomic weight. 
The physical significance of the new variables thus introduced Is 
determined . 


INTRODUCTION 

Shielding is accomplished by interposing a material barrier between 
a source of radiation and the region to be protected. In passing through 
the shield, the particles making up the radiation (which may contain 
particles of zero mass such as photons) undergo collisions with the 
atoms of the shield. In a given collision, a particle may lose some 
portion of its energy while simultaneously changing direction and may 
also generate new particles . The probability of a particular type of 
collision depends upon the energy of the particle before the collision 
and, hence, upon its past history; the number of collisions a particle 
undergoes within the shield depends upon its length of path in the 
shield and, hence, on its original energy and the geometry of its tra- 
jectory. In general, therefore, the particles leaving the shield will 
be of several types and will cover a large range of energies and direc- 
tions . The complexity of the problem makes the solution of the mathe- 
matical equation governing the travel of the particles through the 
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shield (the transport equation) quite difficult. However, the large 
number of possible processes makes the use of composite shields seem 
particularly promising. For example , the first section of a shield 
might slow the particles down sufficiently that the second section could 
absorb them. This means that to find the most efficient shield for a 
given source many co mbi na tions of materials must be considered. Because 
of the difficulties of calculation, direct measurement seems the proper 
procedure for comparing these combinations , but to try all the possi- 
bilities would require a large volume of experimentation. 

In order to avoid both the need for elaborate computation and for 
large-scale experimentation, Bobrowsky (reference l) proposed a method 
(herein designated as s emi empirical) for predicting the effectiveness 
of composite shields on the basis of actual measurements made upon a 
relatively small number of homogeneous shields . The basic idea pre- 
sented in reference 1 is to treat the shield the effect of which is to 
be computed an being built up of a number of layers or elements upon 
ea ch of which the required measurements of shielding effectiveness have 
been made. The results of such measurements serve to relate the distri- 
bution in direction, energy, and surface position of the particles leav- 
ing any element to the distribution in direction, energy , and surface 
position of those entering the element. If the individual elements are 
combined to form a composite shield, the same relations as before hold 
between the distribution of particles entering a particular element and 
the distribution of those leaving that element. Except for the front 
face of the first element and the back face of the last element, however, 
the incident fluxes are unknown. Thus, in the case of n elements and 
p species of particles, there are 2np unknown functions describing 
the emergent fluxes, one for each species of particle at each of the 
two faces of each element. Wien these fluxes are connected with one 
another and with the incident flux, there result 2np integral equa- 
tions which are, however, of a relatively simple structure. 

Briefly, the procedure which has been used by Bobrowsky is to 
approximate the integrals in question by appropriate sums, usually 
with the a ddi tio nal assumption that the incident flux is Independent 
of surface position. The intensity is thus no longer expressed in 
terms of a function of direction and energy but instead is expressed 
in terms of a set of n umb ers, each of which is the average value of the 
intensity within the appropriate range of . energy and of angle. (Actu- 
ally, Bobrowsky suppressed the angular dependence by ass umin g all 
particles to travel in a single d ir ection.) Thus, the set of integral 
equations in the unkno wn intensity functions is replaced by a set of 
ordinary algebraic equations in the intensity numbers . The solution 
of these can be accomplished by ordinary methods. In practice, however, 
the computation may be quite tedious because of the large number of 
unknowns involved) in addition, the amount of experimental data required 
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may still "be excessive. This is especially the case if considerable 
accuracy is desired and, consequently, the intervals of angle and 
energy are narrow. 

The present work, which was conducted at the MCA Lewis laboratory, 
is directed toward extending and simplifying the semi empirical method. 

The first step toward accompli shing these objectives is the determina- 
tion of approximate solutions of the transport equation for two special 
conditions: first, that in which ' multiple elastic scattering dominates; 

and second, that in which the role of multiple scattering with change 
in direction is small compared with other effects. The results obtained 
serve three purposes: (l) They demonstrate the importance of taking 

into account the direction of motion of the particles; (2) they provide 
a means for calculating the effect of any element which is not too thick 
and in which multiple elastic scattering is not too important an effect, 
from the results of measurements on an infinitesimal element of the same 
material; and (3) they suggest a way of simplifying the computations 
necessary in the application of the semiempirical method. 

In order to generalize this method to include the effects of angular 
dependence and to achieve the simplification mentioned, the integral 
equations displaying the angular dependence of the intensity are first 
written down. These equations are to be treated by approximating the 
integrals by sums and solving the resulting algebraic equations. This 
replacement may be accomplished by writing each unknown intensity func- 
tion as the sum of a series of terms, each term being the product of an 
unknown constant coefficient and a known function of the appropriate 
variables. If this representation is substituted into the original 
equations, the integrations may be performed at once, and a set of equa- 
tions in the unknown coefficients results. Knowledge of the approximate 
solutions of the transport equation makes it possible to choose the 
known functions in such a way that all but a few of the terms in each of 
the equations so obtained are small; the resulting equations are there- 
fore amenable to solution by successive approximation. After the coef- 
ficients are obtained, the intensities may be found by summing the 
appropriate series. It is shown, however, that the coefficients them- 
selves contain information of considerable physical significance; for 
some purposes, it may be sufficient to know the coefficients in order 
to estimate the usefulness of the shield. 


ANALYSIS 

Assumptions 

In each case it is assumed that the following conditions hold: 

(l) The elements are homogeneous slabs of finite thickness with 
infinite plane parallel surfaces . 
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(2) No new particles are produced "by collisions. 

(3) The system is in a steady state. 

(4) The inci dent radiation is uniform over the surface and inde- 
pendent of the angle measured around the normal "but not necessarily 
independent of the angle it makes with this normal. 

(5) All sources are external to the elements. 

Conditions (l) and (4) are essential to the type of analysis employed 
here because they represent restrictions on the geometrical aspects of 
the problem -which make the resulting system of equations relatively 
tractable. The requirement with regard to angular distribution contained 
in condition (4) is fulfilled if the sources are disposed in a uniform 
manner with respect to the shield. The assumptions contained in con- 
ditions (2), (3) , and (5) are not essential but are introduced primarily 
to avoid complications in the anal ysis which would obscure the ideas 
presented herein. In particular, reference 1 deals with the case in 
which condition (2) does not hold and reference 2 deals with the case 
in which conditions (2), (3), and (5) do not hold. 


The Individual Element 

The procedure to be followed in simplifying the application of the 
semiempirical method is to determine the way in which the particles 
traverse a single element and to deduce from the results obtained an 
appropriate representation for use in this method. More explicitly, 
the representation sought is in each case to be obtained from an approxi- 
mate solution of the transport equation valid under certain general 
assumptions about the nature of the element. The two cases to be con- 
sidered are that in which scattering dominates and an expansion in terms 
of Legendre polynomials is appropriate (approach l) and that in which 
absorption dom i nates and iteration is appropriate (approach 2). It is 
to be emphasized that the expansions thus obtained are applicable in 
general but are most convenient under the circumstances stated. 

Approach 1: Expansion in Legendre polynomials . - In this section 

the transport equation is transformed by expressing the dependenqe on 
energy in terms of a new variable T] by a Laplace transformation and 
by representing the angular dependence in terms of an expansion in 
Legendre polynomials . It will be shown that for elastic scattering 
the introduction of Tj simplifies the equations considerably, that 
the use of Legendre polynomials is well adapted to the boundary con- 
ditions, and that in the case of multiple scattering only a few terms 
of the exp ans i on need be retained. The need for including angular 
dependence also becomes obvious . 
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The starting point of this discussion is the form of the transport 
equation given by Marshak (reference 3, equation (63); see also appen- 
dix B of this report) here specialized to the time-independent and 
surface -position-independent case with no internal sources: 


X(u)p + i|r(u,p,z) = 

oz 


nu 


du' 


10 


dS’^u^p^z) f([i 0 ^u-u’) h(u’) 


( 1 ) 


All symbols are defined in appendix A. 

Because i{r does not depend upon cp, the right side of equation (l) 
can be integrated over <p ' immediately so that it becomes 



m 


dp , u-u * Mu’ ,p ’, z)h(u* ) 


( 2 ) 


where 

p2* 

fl([i ; |A',u-u') = / f(p 0 ,u-u 1 ) d<p* (3) 

Equation (l) may be further simplified by taking its Laplace trans- 
form with respect to u. If, for purposes of simplicity, it is assumed 
that X(u) and h(u) are constant, then it is possible to put X(u)=l 
(this amounts to a change in the choice of units for X constant) and 
h(u)=l (this amounts to a change in the way f is normalized for con- 
stant h); with these changes the equation becomes (see reference 3, 
equation (17)) 


V- 


+ $ (z,T i,p) = I dp'F(p,p',Tj) § (z,q,p') (4) 


-1 


where $ and F are defined by 


= I cLu e -T1U t(u,p,z) 
0 


( 5 ) 


and 
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F(p,p*,Tl) = du e" 7 ^ fq(p,p',u ') 
0 


( 6 ) 


respectively . In equation (4) the integration over u' has "been per- 
formed hy making use of the fact that u’ enters fq OD dy through 
u-u‘ . The extension to the case -where X and h* are not constant is 
imm ediate (reference 3, equation (80)). It Is only necessary to expand 
each of these quantities in a po-wer series in E: 



(?) 

( 8 ) 


-where the quantities Aj and Hj are constants and the j 1 s represent 

some set of numbers of either sign which may be, but are not necessarily, 
integers. These expansions need contain only a few terms if X and 
h do not vary too fast with u. The representation given is not very 
adequate if resonance is marked; but the error is not serious for a 
narrow resonance band because the chance of a particle attaining an 
energy within such a band is small. If equations (7) and (8) are sub- 
stituted in equation (l) and the Laplace transform is obtained, equa- 
tion (l) becomes 


r 





^(z,T]+j,p.) 

dz 


+ * (z,Tbn) = 


I 


pi 



dp.’ 


F(p,pSil) $ (z,Tj+j,p») 


0) 


It is to be emphasized that the validity of equation (9) depends 
upon the assumption that u and u* enter fq only through the com- 
bination u-u* . This restriction on fq, while not necessarily always 
holding, is correct in many cases of physical significance (reference 3, 
p. 186), for example, whenever the scattering can be represented as a 
classical collision between the incident particle and a free particle 
at rest. The same result holds in quantum mechanics for the case of 
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so-called S-wave scattering, tliat is, the case when the "wavelength" of 
the incident particle is great compared with the range of the scattering 
force. In the general case, the function fp can he written in the 
form 

f l = f i( u " u S u ') ( 10 ) 

so that the treatment of this case can he extended hy expanding fp in 

a series in e just as was done in the case where • h(u) was not 

constant . 

In order to couplet e the formulation of the problem-, it is neces- 
sary to specify the boundary conditions. This may he done quite simply 
as follows . Let z = 0 represent the face upon which the radiation is 
incident and let z = Z represent the opposite face. Then the condi- 
tions in question are 

= A(tj,p) p > 0 

z=0 


. =0 P < 0 

I z=Z 

The quantity A(t],p) is here an arbitrary known function of the 
appropriate variables. Physically, these equations mean that the flux 
incident upon the front face of the element is known while the flux 
incident upon the back face is zero. An. individual element of the shield 
is shown in the following diagram: 


$(Z>T),p) = 0 
for p < 0 

z=Z 




z=0 
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Equations (4) as they stand are still integral equations with, respect 
to p as well as differential equations with respect to z. It is 
desirable to transform from this set of integro -differential equations 
to a (infinite) set of ordinary linear differential equations. This 
can he simply done if and F are expanded in terms of complete set 
of orthogonal functions of p. Marshal; (reference 3, equation (65a)) 
used for this purpose the set of Legendre polynomials of p. For the 
present problem this is inconvenient because the boundary conditions 
(equation (ll)) impose a sharp distinction between forward (p > 0) 
and backward (p < o) directions. It is therefore natural to make use 
rather of an expansion in which this dichotomy is reflected in the 
expansion of the function § in terms of the variable p. A procedure 
in which this is accomplished and in which some of the advantages of 
the use of the Legendre polynomials are retained is based on the use 
of the representation: 


$(z,il,p) 



(2n+l)<p£0i,z)P n (2p-l) 


( 2n+l) (f^ (q , z ) P ( -2p-l) 


p > 0 


p < 0 


( 12 ) 


The + and - superscripts on the expansion coefficients <p 
indicate, of course, that they are associated, respectively, with posi- 
tive and negative values of p and hence with forward and backward 
directions . The coefficients have thus a significance different from 
that to be associated with the usual expansion in terms of Legendre 
poly nom ial s (reference 4) . It is to be emphasized, however, that the 
number of nodes of the Legendre polynomial of degree n increases as 
n increases; and, hence; the larger is n, the greater is the degree 
of angular assymetry to be associated with q> n . It will be advantageous 
for future work to define 


P n (2p-1) =< 


P n (2p-1) 0 < p < 1 


lo 


otherwise 


and 


’^nC -2 ^ -1 ) -1 < P < 0 


P’(-2p-l) =< 


.0 


otherwise 


( 13 ) 
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With this definition § may he -written 


§ 


= £ |( 2tl+X) | jpJ(T] ,, z ) P^C 2|n-l) + <p"(ri,z)p~(-2q-l)^ (14) 


It is to he noted that two P*s for -which the n's and the super- 
scripts are not hoth identical are orthogonal. This is one of the chief 
advantages of the expansion. Because of the isotropy of the element , it 
is assumed that f depends on 6 and 0' only through |pq| and, 
therefore, F may he -written in' terms of these functions as 


F = 


£,(2**1) [pJ(2ia-l)P+ t (2^’ -1) + P~(-2n-l)p",(-2^-l)J 
f+^ t [p+(2p-l)p-,(-2n«-l) + p-(-2|a-l)P+ t (2n«-l)j| 


(15) 


This equation serves to define the expansion coefficients f . (These 
coefficients have the property that f n n , = f nI n .) When equations (14) 
and (15) are substituted into the right side of equation (4) , there results 


Si . - 


= E ( 2n+1 ) j J'2,r<,P+(2n-l) + T n ,p-(-2|i-l)l 

n,n’ L L. J 


(16) 


where use is' made again of the properties of the Legendre polynomials 
(see reference 4) . If equation (16) is multiplied hy PE(±2q-l) and 
the resulting equation integrated over the range -1 < q < 1, the fol- 
lowing set of differential equations is obtained: 


1 <^0 1 ^1 ± Y"T + + ± _+ - J*] 

2 Sz ± 2 + ^0 “ £/|5o, m % + f 0,m ^mj 


, 1 K x 1 n+1 K+l x 1 a a ^-l , * 
2 bz 2 2n+l Sz 2 2n+l Sz n 




Yfr cp* + cp^l n > 0 

£nm m nm mj 


( 17 ) 


J 
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To tills set are to "be appended the boundary conditions -which here take 
the simple form 


^n z=0 




= 0 


n |z=Z 

-where Z is the thickness of the element and 

ni 


A a (Tl) = 


0 


p*( 2 p-i) $ ( z ^^) 


(18) 


2=0 


dp 


From the viewpoint of shielding, the problem is to find | and 

^n|z=Z ^ ' berms of tlle A n( T l) * Generally speaking, for large Z one • 

may expect these quantities to decrease as n increases in the case 
where there is little absorption because multiple scattering enhances 
the symmetry. A simple illustration of this is given by the following 
example: 


Example of approach 1: Scattering of neutrons by a thick "heavy 

shield" without absorption . - When a neutron is scattered by an element 
of relatively high atomic weight, the neutron loses very little energy 
and the angular distribution of scattered neutrons is nearly uniform. 

As an approximation, the change in energy may be taken as zero and the 
angular distribution as uniform so that F is constant. As it is 
assumed that there is no absorption, F is normal ized in such a way 
that 



(19) 


(see equations (3) and (6) and equation (Al) in appendix A) . Thus it 

follows that, to the approximation taken here, F = i. With the con- 

2 


sistent neglect of all <P n *s for -which n exceeds 1, the differential 
equations (17) for become 


1 ^1 1 S<P 0 
± 2 S z ± 2 d z 


± 1 ± 1 T 

+ ^0 = 2 + 2 ^0 


1 *$ 

± v ^ — ± 


&P' 


3 dz Sz 


1 + 2cpJ =0 


. ( 20 ) 
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The boundary conditions are 

> ( 21 ) 

J 

The dependence on T] is suppressed because the value of <P n for a 
particular T) is not related to its value for any other T). These 
equations , being a set of linear differential equations with constant 
coef f icients , may be solved by the usual methods. As remarked pre- 
viously, the feature of interest here is the way in which the element 
imposes angular uniformity upon the incident radiation when Z is large 
compared with 1 (that is , with X) . This effect is demonstrated in the 
following approximation to the solution in which all but leading terms 
are omitted. With such an approximation, the result may be written as 



*8)2=0 = ^o - ( 6 -W 3 ) \ 

*l|z=0 = (7-W3) A x 


( 22 ) 


The omitted terms are of the order of A^ (or Ap) multiplied by l/Z. 
Because $ must be positive at z=0, it follows from equation (l4) 
that Aj_ must be less than one-third Aq. If it is noted that 7-4'\/3 
is about 0.1 and is about 0.9, it follows that 

is considerably smaller than To) 2 = 0 * Because q>Q represents a uni form 
emergent flux and (p£ is associated with radiation displaying a certain 
degree of asymmetry, this confirms the assertion that the shield acts 
to promote the uniformity of the radiation. If the next stage of the 
calculation performed is employed by assuming that cp^ = 0 for n > 3 
and substituting the results of the present approximation into the 
equations thus obtained for <Pg and <Pg, it follows that, in accord 
with the original, appro xim ation, these last are also small. 


The requirement that there be no absorption is reflected in the 
value of F. In particular, if the absorption were not zero, F would 
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"be less than l/2 and <pq would he decreased in magnitude relative to 
the asymmetric terms -whose importance would therefore increase. In 
addition, for an element the constituents of which are of lower atomic 
weight, the coefficients f m i do not vanish for n' greater than 
zero and they also vary with tj . The effect of this dependence is to 
increase the number of scatterings required to produce the same amount 
of uniformity. 

Under the heading of approach 1, a procedure has been developed for 
dealing with an element in which scattering is the do minant effect. In 
particular, it has been assumed that X and h are constant, or nearly 
so, and that the scattering process is of the S-wave type. Under these 
circumstances it follows thjt an appropriate representation is in terms 
of the intensity numbers <P n (T]) . The expansion thus obtained takes 
account of the sharp distinction between the forward and backward direc- 
tions, gives explicit recognition to the expected effects of scattering 
in producing angular uniformity by segregating terms implying different 
degrees of angular asymmetry, and takes advantage of the simplification 
resulting from the use of the Laplace transformation. 

Approach 2: Iteration . - Although the expansion discussed in 

approach 1 may be used in any problem, its practical application is 
primarily found in those cases where there is elastic scattering and 
where the element imposes a certain degree of angular uniformity upon 
the emergent radiation. In the contrary case where the emergent radia- 
tion retains essentially the nonuniformity of the incident flux or 
where such nonuniformity is imposed, it is convenient to use a formalism 
in which the intensity is expressed directly as a function of the angle 
and the energy. 

This will be the case if multiple scattering with change in direc- 
tion is not an important process. The procedure used here is an iter- 
ative one which again shows the importance of including angular depend- 
ence and which permits direct solution of the transport equation for a 
single element under the circumstances indicated. 

The relevant transport equation with no scattering leading to 
directional changes is (see appendix B) 

p + Q(E) llr(E) =|q 1 (E,E I ) t(E') dE’ (23) 

Here Q is l/k and q-j_(E,E’ ) dE is the proportion of radiation of 
energy E* degraded into the range between E and E+dE. The func- 
tion q^_ vanishes for E > E 1 (no gain of energy) . It may be noted 
that, as this equation stands, p enters only as a parameter. If i{r 
is written as 
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S3 

CD 

CD 


^oo 

i}r(E,ii,z) = T(E,z,E') ^(ES^O) dE* (24) 

uo 


egu.ati.on (23) is satisfied if T fulfills the boundary condition 
T(E,Q,E’) = 5(E-E') and the equation 


+ q(e) t(E,z,E') = q-i (E,E") T(E", z,E> 
oz 


) dE" (25) 


E' and p are parameters in this equation. 

An approximate solution of equation (25) may now he obtained by 
the use of iteration. This approach is especially appropriate because 
it includes as a special case the method of breaking the energy range 
into intervals and replaning the integral on the right side by a sum 
essentially employed in reference 1. In order to carry ’out the itera- 
tion, T is expressed in the form 

oo 

T = ^ (26) 

j=0 

where T^) satisfies the boundary condition 

O 

T^j) (E,0,E*) = &JQ 5(E-E') (27) 


and represents the contribution to T of particles that have started 
out with energy E* and have been scattered j times. It follows that 

T^) is of the order of q^ because the latter is essentially the 
probability of j scatterings. Consequently, if equation (26) is sub- 
stituted into equation (27) and terms of the same order in q-^ are equated, 
there results 


p — — + Q(e) t(°)(E,z,E') = 0 

oz 


-\ 


,(j) 


r®max 


— { E ; Z ; E . ') + Q(E) T^(E,z,E’) = 


u 


0 


(E,E") T^ _1 )(E",z,E’) dE" 

(28) 


& 
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the solution of -which, subject to the boundary condition (27), is 
(appendix C) 

*\ 


t(°)(E,z,E') = 5 (E-E * ) 




(E,z,E«) 


^(E'E') f -Q(E«)^- -Q(E]~ 

— T— r — r- e R -e M- 

Q(E)-Q(E') l 


T (j) (e^°^z,E^^ 


) 



^(E^V^X 


) (29) 



e-Q(E^ e ^)z/p j " 1 
Q(E^ r ^)-Q(E^ e ^) k=l 



j > 1 


J 


In case tiro of the Q's are equal, the expression is to be evaluated 
by initially treating them as distinct, grouping all terms in -which 
their difference is contained in the denom in ator, and then taking the 
li mit. In case the set of energies which the particle can lose upon 
scattering is discrete or in case the integral is approximated by a 
sum, the function q^ can be represented as a sum of 6 functions. 
If this sum contains t terms, then the series for T breaks off 
after T^) . (This depends upon the fact that q^(E,E’) is zero for 
E larger than E * . ) 


A variation of the preceding treatment can be applied to the case 
■where there exists a scattering process which produces a change in both 
direction and energy as well as the original process which affects only 
the energy. In this case the differential equation may be written 
(appendix B) 


H + e(E) ♦(*,„,*) = 


'^max 


0 


q 1 (E,E I ) *(E*,p,z) dE* + 


PI 


U-i u 


"®max 


cl 2 (E,p;E',h«) *(E*,p*,z) dE* dp* (30) 
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where qg^pjE'^p’) dE ap is the proportion of particles with energy E* 
p.nd direction p' scattered as a result of the second process so as to 
have energies in the range between E and E+dE and directions in the 
range between p and p+dp. Just as was done in equation (24) i{r may 
be expressed in the form 


*(E,P,z) 



[^(E^p^ZjE'jpO^CES^S 0 ) clE’ <3p* 


where T satisfies the boundary conditions 


(31) 


T(E,p,0;E l ,p') = B(p-p’) 5(E-E») p > 0 
T(E,p,Z;E* ,p’) =0 p < 0 


(32) 


and the equation 


^ ^T(E,p,zjE»,pQ 
dz 


+ Q(E) T(E,p,z;E*,p*) = 


q 1 (E,E") T(E",p,z;E',p*) dE” + q 2 (E,p,E",p") T(E",p",z$E« ,p* ) dE" dp” 

J w (33) 

If the second integral on the right side of equation (33) is small 
in comparison with either the second term on the left side of equa- 
tion (33) or the first integral on the right side of equation (33) , then 
the approximate solution (equation (29)) obtained for equation (25) may 
be adapted for use as the first approximation for T by applying the 
boundary conditions (32) . The resulting function, which may be denoted 

as T^^EjPjZjE'jP’),, satisfies 

p — + q(e) t(°)(e,p,z$e»,p0 = 

oz 



(E,E") T^ 


(E",p,z;E 


,p') dE" 


and the boundary conditions 

T( 0 )(E",p,0jE , ,p’) = B(p-p') 8 (E-E * ) 

T^°^(E",p,Z)E',p) =0 p <0 • 


( 34 ) 


p > 0 
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A "better approximation is obtained "by putting 

T = + T (■*-) 

where is as "before and 


(35) 


satisfies the conditions 


T^CEjPjOjEJp’) =0 p> 0 

T^CE^ZjE’ =0 p< 0 

If equation ( 35 ) is substituted into equation (33) and if the product 

-(l) -(0) 

q.gT i6 neglected in comparison with q_gl , the following equation 

is obtained for : 

P — — - Z ; E - ^ * )- + Q(E) T^^PjZjESp’) = 
dz 


q 1 (E,E") (E",p,z$E’ ,p‘) dE" + 


<j^(E,p$E",p n ) T^(E",(i' , ; z;E , ; |i’) dE" dp" 


(36) 


This equation in turn may be treated by iteration. If the symbol 


S(E,p,z;E‘,P*) = 


q 2 (E,p E",p") (E",p",zjE , p I ) dE" dp" 


is introduced, a first approximation to T^, obtained by neglecting 
the first integral on the right side, is given by 


T^ 1 A \e,p,z;E’,p>) = 




r\ z 


^ Jo 


exp 


Q(s)(E-z) 


S(E,p,£;E«,p«) d£ for p > 0 


and 


r>Z’ 


1 


exp 


<»(*0(£-g) 


P -* 


S(E,p,£;E‘,p«) a£ for p < 0 


( 37 ) 


J 
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This may then he substituted into the first integral on the right side 

of equation (36) and a better approximation for T^ obtained. The 
iteration may, of course, be continued along these general lines, but 
its usefulness is, in general, greatly diminished as the length of the 
process is increased. In any case, the early appropriation •will be good 
if the contribution from the integral involving qg is not too large; 
that is, if multiple scattering with angular deviation does not occur 
to any great extent. 

It is apparent that in the two cases considered (q^ small or 
essentially zero) the emergent flux will depend strongly on angle for 
an element of reasonable thickness; indeed, for a sufficiently thick 
absorber, it follows from the form of that this flux is mainly 

concentrated within a cone of half -angle of order \J2.JQL where the 
value of Q to be used is a typical one. Consequently, an analysis 
in which the angular dependence is given by specifying the amount of 
radiation within each of a set of ranges into whi ch p is divided is 
appropriate . 

The results obtained here serve to elucidate the case in which 
scattering with angular deviation is not a very significant process so 
that it may be treated as no more than a perturbation. An exact solu- 
tion is obtained by the iteration method for the case where the direc- 
tion of motion of the particles is unchanged; but, unless Q and q^ 
have a particularly simple form, its application is obviously inconvenient 
where there is much multiple scattering for the reason that, under such 
circumstances, a large number of terms must be retained in the series 
for T. The process for dealing with the perturbing effects of scatter- 
ing with angular deviation is described. 


Semiempirical Method 

The results of the previous section, entitled "Individual Element, 11 
find application to the treatment of shields by a modification of the 
semiempirical method. In this formalism, a shield intercepting a flux 
of particles is to be represented in the following fashion. Radiation 
of an intensity independent of position on the face falls upon the first 
of a set of n plane parallel elements. As a result, there is incident 
upon both faces of each element a flux of radiation; consequently, there 
emerges from each face a certain quantity of radiation. Let J^(E,p) 

represent the number of particles with energy E and direction p per 
unit range of E per unit range of p per unit time incident upon each 

unit area of the i^ element. Similarly, let I^(E,p) represent the 
flux of radiation emergent from the i^ element. The following dia- 

_l_T_ 

gram depicts the l element of the shield: 
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Positive p indicates radiation incident upon the front face, and neg- 
ative p, that incident upon the rear face. For I, then, negative p 
is associated -with radiation leaving the front face and positive p, 
that leaving the rear face. The manner in -which the presence of the 
element modifies the radiation is given hy the equation 


Ii(E,p) = 


E^(E,p;p i ,E^) J ± (p i ,E i ) <% dp ± 


(38) 


where K i (E,pjPjL,E^) dE dp represents the proportion of particles with 
energy E^ and direction p^ which is modified hy the i ^ 1 element 
so as to emerge with energy between E and E+dE and direction between 
p and p+dp . 


In addition, the I’s and J’s satisfy the following equations: 

'N 


I ± (E,p) = J ±+1 (E,p) p > 0 

J ± (E,p) = I i+1 (E,p) p < 0 


J n (E,p) =0 p < 0 


p < 0 


> ( 39 ) 
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and I-^(E,p) is known for p > 0. It follows from symmetry considera- 
tions that 


KiCEjPjE^Pi) = Kj^Ej-pjE^-Pi) 


(40) 


The previous set of equations Is, as noted in the introduction, best 
treated by some method of approximation. In reference 1 this was done 
by replacing the integral by a sum; as shown previously (see remarks 
following equation (37)), this is appropriate where the dominant proc- 
esses are absorption and loss of energy without scattering. 


However, in the case when elastic scattering is of primary import- 
ance, another procedure suggests itself (see Approach l) . In this case 
it is natural to expand I and J in terms of the polynomials 
P n (.±2p-l) and to use the Laplace transforms of the coefficients so 


obtained as the intensity numbers . This may be done by defining the 
quantities 


no 


*£n(0 « 


ol 


®max'J 


\ 


Pj(2p-1) e -T * U I ± (E,p) dp dE 


0 


r0 


iIOi) = 


no 


/ (41) 


®maxU 


P*(-2p-l) e -T1U I ± (E,p) dp dE 

-1 J 


with analogous expressions for and J“. In terms of these quanti- 

ties, equations (38) and (39) become, respectively, when \ is inde- 
pendent of E and there is no absorption. 


and 



- -’Ll, 1 * n 
1 ■= 11 

j;0i) = o 


*> 

> 

J 


( 43 ) 



/ 
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■where the quantities k !^ t (tj) and k l|~~n t (~n) are defined hy: 

^iSn’C^) = 


(2n*+l) 



n 


oj 


pJ(2h-1)P+,(2ki ± - 1) e~^ U K, (E^E^m ) dp du* dE 

(44) 


and 


^wOl) = 


(2n‘+l) 


po 

'XL 

®max'° 

Ov. 


00 


pJ(2h-1)P-,(-2h 1 -1) e‘ nU K^E^E^^) dp dE 


The equations here contain tj only as a parameter; hecauBe the element 
has the effect of producing a distribution of intensity nearly inde- 
pendent of p, it should he sufficient to^eonsider only a few values 
of n. In case the mean free path is not constant, the right side will 
depend on other values of t) besides the one appearing on the left 
side (see equation (7)). In this case, if the average rate of varia- 
tion over a wide region is not too large, it is possible to apply suc- 
cessive approximation by first omitting these extra terms . If this is 
impractical, the best procedure is simply to expand I and J in 
terms of Legendre polynomials as before, but with the coefficients func- 
tions of energy rather than of tj. The integrals may then be replaced 
by sums as in reference 1. 

Even in the case in ■which the expression in terms of the Laplace 
transform expedites the solution, there remains the problem of Inter- 
preting the results obtained, which will give quantities like (D) • 

The corresponding function of u and p may be obt ain ed by performing 
the inverse Laplace transformation with regard to tj, multiplying the 
resulting function of u by (2n+l)pJ, and s umming over n. Each of 
these procedures, and particularly the first, is likely to be quite • 
tedious. It is therefore desirable to avoid the 'actual calculation; 
and this may sometimes_be done because, as the following considerations 
show, the quantities I^ n ( T j) themselves contain considerable informa- 
tion of physical sig nif icance. If the dependence on n and p Is sup- 
pressed, the relation between I and I is given by 


9P£Z 
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I(f]) 


= J e _T * U I(u) du 


(45) 


Here l(u) may be taken as representing the rate of flow of particles 
per uni t area per unit time per unit range of u. In view of the con- 
nection between u and E, I ( rj ) may be rewritten in the form 



(46) 


Thus., l(o)_ is the rate of flow of particles per unit area per 
unit time, E max l(l) is the average rate of flow of energy per unit area 
per unit time, etc. The quantity l(i}) therefore finds interpretation 
in terms of particle current density, of energy current density, and of 
the moments of the energy flow. If l(u) itself must be known, the 
values of u of interest will usually be large because the important 
question will be whether the particles have been sufficiently slowed 
down to reach an absorption band in a succeeding element. In this case, 
use may be made of known methods for approximating the inverse trans- 
formation for large u (see reference 3, p. 204). If the dependence 
on n and p Is now considered, it follows in a similar fashion that 
much may be deduced about the directional flow of particles directly 
from I n . Thus the total current in the forward direction is (from 
equation (46)): 


nO 


^^max 


o 1 m 

■ L 


l(E,p)p dp dE 
0 n 


(2n+l) ij(o) 

V 

= I |jo(°> + 


0 


P n (2p-l)p dp 


(47) 


It Is therefore possible to learn a good deal about the particle distri- 
bution without transforming from the variables (n,7]) to the vari- 
ables (p,E) . 


Up to now, the results of the section "Individual Element" have 
been considered from the viewpoint of their significance in the solution 
of a problem involving a number of elements whose separate effects upon 
the incident radiation have been determined experimentally. However, 
the results obtained in approach 2 of that section may also be used to 
determine, to a good approximation, the effect of a thin element of a 
given material for which the constants associated with each possible 
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process are known. The semi empirical method may then he used to obtain 
the required data for thick shields. Because the term "thin" , as used 
herein, means only sufficiently thin that multiple scattering need not 
he considered, it -would appear that this procedure -would provide a con- 
venient and reasonably flexible -way of obtaining the needed constants 
for a variety of elements. 


CONCLUDING REMARKS 

The problem of a layered shield consisting of a number of plane 
parallel elements of infinite extent but arbitrary thickness has been 
treated. By use of the transport equation, an explicit approximate 
solution has been obtained for a single element for the case -where mul- 
tiple scattering with angular deviation is not dominant and the form 
of the solution investigated in the case where multiple scattering with 
angular deviation is the principal effect . In each case, an expansion 
is given for the angular dependence ' which is well suited to the boundary 
conditions of the problem and -is rapidly convergent . It is pointed out 
that the use of an expansion of this sort and of the Laplace transform 
with regard to energy (where scattering with angular deviation is dom- 
inant) should considerably reduce the labor in a semiempirical treat- 
ment of shielding. On the other hand, the semiempirical approach may 
be used to extend the single scattering result to a thicker element 
where multiple scattering is of importance. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 8, 1950 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: 


AnC^) 

A(t],p) 

E 

®max 

F(p,p’,ri) 

1 

f(MQ,u-u*) 


f^.Cri) 

nn T v 1 


H(u) 

li(E^) 


expansion coefficient for A(t],p) in series 

co 

A(t],p) (2n+l)A n (t])P n (2p.-l) 

n=0 

value of $ at z = 0 for p > 0 
energy of a particle 

maximum energy of any particle entering shield 
Laplace transform with respect to u of f-^p^p* ,u) 


r^2si 


0 


f(p Q ,u-u'} d<p* 


proportion of particles with direction and 

energy E' scattered through the angle 0 into the 
direction (0^<p) with energy E per unit solid angle 
per unit range of u 

expansion coefficient for expression of F in terms 
of P+(2p-l), P+,(2p-l), P~(-2p-l), and P" f (-2p-l) 
divided hy ( 2n+l) 

coefficient of e in expansion h(u) = YY e 




•MeL 

K( u ) 


<3 


coefficient of E^ in expansion h(u) = ]T^ kjE 


number of particles per unit time per unit area leaving 


ith 


element of shield with energy E and angle 


cos - ''' p per unit range of energy and unit range of p 
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4aOl) 


J ± (E,n) 




OO 


m 


^nifly v 

no 

^iis.y < 


0 

no 


? n (2\i-l)e~ W I-jfop) dp dE 


P n (-2p-l)e _11la I ± (E,p) dp. dE 


number of particles per unit time per unit area incident 
upon i^ element of shield with energy E and 
angle cos - ^ p per unit range of energy and unit 
range of p 


U 


^ 0 pi 

®maxv- ® 


PnC 2 ^- 1 ) 


-■Hp 

e 


J ± (E,p) 


dp dE 


no 


no 


JiXL^) 


-^niaxvJ ' 


P n (-2p-l)e" TW1 J ± (E,p) dp dE 


K^E^p 1 jp^E^) proportion of particles incident upon i^ 1 element 

whose original energy is E^ and direction, p^ 
scattered into energy E per unit range of E and 
direction p per unit range of p 


% nn'M 




no 


= (2n4-l) 


= (2n'+l) 


m 




m 


r %0 

nl 


o^ 


oUo 
ro 

-l 


F n (2p-l)P n ,(2p-l)e- , l u <% dp dE 


P n (2p-l)P n ,(-2p 1 -l)e _,ltl g,(E,p ? B^„,p,) dp* dp dE 


P n ( x ) 


Legendre polynomial of degree n 
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p*(2h-i) 

Q(E) 

q^E') 

q 2 (E,p,E»,p‘) 

S(Ej(J-^z;E’ ,p’) 

T(E,z,E*) 

T(E,ii,zjE , ,ii») 

t( j) 

5p(j) 

u 

Z 

z 

&(x) 


' P n (2n-1) Ji > 0 

p < 0 



(-2p-l) 


p < 0 


p > 0 


X( u ) 


proportion of particles -with energy E* scattered into 
state with energy E per uni t range of E 

proportion of particles with angle cos - -*- p* and 

energy E 1 scattered into state with angle cos - -* - p 
and energy E per unit range of E per unit range 
of p 


C PP^max 


q. 2 (E,p;E ,, ,p")T 0 (E%p", Z ;E’,p>) dE" dp" 


J-iU 


o 


proportion of particles with energy E 1 incident upon 
face of element reaching energy E per unit range E 
upon penetrating to depth z for g 2 = 0 

proportion of particles with energy E' and angle 
cos -1 p* incident upon surface of layer reaching 
energy E per unit range of E and angle cos - ^- p 
per unit range of p upon penetrating to depth z 

.P* 1 approximation to T 

"ttl 

j approximation to T 

E 

n max 


thickness of layer 

distance of point inside layer from the front face of 
the layer 

Direc delta function 



26 


HACA TH 2647 


j Kronecker delta 

H variable of Laplace transform 


0 

0 



X s ( u ) 




qfk 

n. 


ifr(E,p,z) 


i|r(u,n,z) 


2 


angle "between direction of scattered particle and direc- 
tion of particle "before scattering 

colatitude measured with respect to normal to face of 
shield upon which radiation is incident 

coefficient of e“^ u in expansion X(u) = Y1 A. e"J u 
total mean free path 

coefficient of E^ in expansion X(u) = X* E^ 

scattering mean free path 
cosine 0 
cosine 0 

Laplace transform of tJt with respect to u 
longitude about normal to surface 


expansion coefficient in expression for $ in terms of 
(2n+l) P$(±2p-l) 

number of particles per unit time per unit area per 
unit range of E per unit solid angle of energy E 
and direction cos -X p and location z divided by 
X(u) 

number of particles per unit time per unit area per 
unit range of u per unit solid, angle of energy E, 
direction cos -X p, and location z divided by X(u) 

solid angle 


All Integrals for which no limits are specified are to be extended 
over the entire range of the variable concerned. 
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APPENDIX B 


DERIVATION AND SIGNIFICANCE OF TRANSPORT EQUATION 

The transport equation simply expresses the fact that any particle 
entering a given region must., in the steady state, he absorbed, be scat- 
tered, or else leave the region unchanged, and that any particle leaving,, 
the region must, in the steady state "with no internal sources, either 
be replaced by a like particle entering the region or else be produced 
■when a particle "with a different momentum is scattered. In order to 
state this in mathematical terms, it is necessary only to consider the 
balance within a region R bounded by a surface S. The following 
contributions and losses occur: 

The number of particles with energies between E and E+dE (where 
E = E mnir e~ u and dE = - E du) movi^ in directions lying within a cone 
of solid angle <32 around the direction specified by the unit vector n 
which enter R through the surface S per unit time is 


\i|r n * dS du dfi 


where cL§ is a vector element of area. This holds because is 

the vector current density of particles per unit time per unit area per 
unit solid angle per unit range of u. 

The rate at which particles of this energy and direction disappear 
within R is 


ijr dv du <32 


where dv is an element of volume. This follows from the fact that 
the linear rate of decrease of the density of particles is simply the 
density divided by the mean free path and that the time rate of decrease 
of density is equal to the linear rate multiplied by the speed; but this 
last product is just ijr. 

The rate at which particles with other momenta are scattered within 
R so as to achieve an energy between E and E+dE and a direction 
within the solid angle <32 centered around n is 
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rvu 


du 1 1 

Uo w 


d£2’ 


dv Tjr(u’ ^tJ-SzifCp-Q^U-uOtLCu.’) ClU 


This may "be seen as follows: h(u’) f(|OQ,u-u’) du <12 is the proportion 

of particles either absorbed or scattered which originally have an energy 
*E' and a direction n* and which are scattered so as to achieve an 
energy lying in the range between E and E4dE and a direction lying 
within the solid angle dS2 centered around n. Thus when thiB product 
is multiplied by i{r and integrated over a ll directions dS2’ , over u’ 
between 0 and u (f=0 for u < u' because of the assumed degrada- 
tion of energy), and over R, the desired rate is obtained. A few other 
remarks about f may be in order: (l) Because every scattered par- 

ticle must attain some direction and energy, f satisfies the equation 




u 1 



dS2 f (pq,u-u* ) = 1 


(Bl) 


(2) The dependence of f upon u-u' results from the assumption of 
classical elastic or quantum mechanical S-wave scattering (that is, 
scattering symmetrical with regard to a coordinate system located at 
the center of gravity of the system participating in the scattering 
collision). (See reference 3j pp. 187, 188.) 


The resulting equation for balance is: 



X(u) t(ii,p,z) n 




^^(u'^Sz) f(|x 0 ,u-u>) h(u’) 


(B2) 


If the left side of equation (B2) is converted to a volume integral by 
Gauss's theorem and the integrands in the integrals over R occurring 
in the equation set equal in view of the fact that R is arbitrary, 
the following result is obtained: 

X(u) n • Vt(u,p,z) + t 


^'^(u-Sm-Sz) f(n 0 ,u-u I ) h(u’) (B3) 


du’ 


0 
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Because ijr depends on position only through z and because the com- 
ponent of n in the z direction is cos 6, n^v^ reduces to 

cos 9 ^ and equation (B3) becomes identical -with equation (l) . 
oz 


Similarly, equation (23) is obtained by writing ijr as a function 
of E rather than u and by assuming that a particle does not change 


f(po,u-u') h(u 


’) du' becomes 


direction on scattering so that 

n 

q 1 (E,E*)S(p-p‘) dE* 

. Equation (20) is the same as equation (23) 


2 it 


except that in equation (20) Jf (pq,u-u* ) h(u’) du* is replaced by 

the more general function 


r<- 


q 1 (E,E’)5(^ , )+l 2 ( E 


2 it 


jH?E«,n')1 


dE’ 


The q-^ term represents all the scattering without angular deviation, 
as before, and the q 2 term allows for an additional scattering process 
which produaes a change in direction. 
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APPENDIX C 

DERIVATION OF EQUATIONS (29) 

It may be verified "by direct substitution that the expressions for 

T (0) 

and t(i) 

'given hy equations (29) satisfy equation (28) subject 
to condition (27) . For values of j greater than one, the correctness 
of equations (29) may he checked hy induction. Because the general 

expression given for T^ is correct for j =1, it remains onlv to 
show that, if equations (27) and (28) are used to compute T^ + ^-' on 
the assumption that T^^ is given hy equations (29), then the result- 
ing expression for T^ agrees with that obtained from equa- 

tions (29). The formal solution of equations (27) and (28) for t(J + ^) 
with given is 


'^E nrax ryz. 

Jo Jo 


c -|ci(e(J +1 ))z - q(e(J +1 ))5] I 


t(J)(e(°),£,eM} at; dE^ 

(Cl) 


Here T^ stands for the function specified in equations (29); when the 
explicit substitution is performed, the right side of equation (Cl) 
takes the form of a multiple integral with respect to dEC 1 ) . . . dE(«5 -1 ) . 
The integrand of this multiple integral may he written, upon integration 
with respect to £, as 


Fornax 


» s=0 1 n=0 ‘r=0 q7e^ 


-q(eCb))1 


i^n 


Q(E^) -Q(E^) 


+ £ n 

n=0 r=0 
rj&. 


r Q( E (j+l))^ 


Q(E^)- Q(E ( - n ' > ) Q(E^) - Q(E^ +1 ^) 


(C2) 


In order that this result agree with equation (29) it is necessary 
only that the following relation hold: 
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<1 

E n 

n=0 r= : 


Q(E^ r ^) -Q(E^ n ^) Q(E^) - Q(E^ +1 ^)j 


1=0 Q(E^) - Q(E^ +1 ^) 


( 03 ) 


Equation (C3), in turn. 


is valid if the following expression van i shes: 


£ £ Q(E( r ))-Q(E(3 +1 )) _ 1 

n=° r=0 Q(E^ r ^) - Q(E^) 

rfn 


(04) 


Equation (C4), however, is just a polynomial of degree (j-l) or less 
in Q(e(< 5 +1 )) which vanishes whenever q(e(< 5 +1 )) takes on one of the 

values q(e(°)), . . . , q(e( 3)) and, hence must he identically zero. 
This completes the verification of equations (29) . 
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